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Abstract 

Discerning the expiration status (non-rejected and rejected) of edible vegetable oils is very significant because of the 

hazardous primary and secondary oxidation products. Therefore, it is of outmost importance to monitor the quality and 
safety of these oils. Based on previous literature, reports and experimental observation, the oil color changes during 
oxidation. Thus, the present study investigates the use of image processing and linear discriminant analysis (LDA) for the 
classification of non-rejected and rejected edible vegetable oils during oxidation process at 85°C, with respect to the 
induced period in both primary and secondary oxidation of four oil type (Olive, Sunflower, Palm and Soybean). The 
purpose of this study was to find less costly and quicker methods with environmental protection, by using the color spaces 
(RGB, HSI, L*a*b* with Grayscale) instead of chemical analyses to determine the expiration status of edible vegetable 
oils. Results of this study indicated that the best classification for expiration status of known oils according to induced 
period of peroxide value in each color space, was achieved with LDA model were for palm with 100% (HSI and 
Grayscale), olive with 84.61% (L*a*b* and RGB), soybean with 95% (Grayscale) and sunflower with 100% (RGB and 
HSD), also in induced period of carbonyl value test, the best classification performance was achieved in palm with 100% 
(L*a*b*), olive with 100% (L*a*b*), soybean with 89.47% and sunflower with 95% (HSI). 


Keywords: Edible Vegetable Oil; Oxidation; Peroxide Value; Carbonyl Value; Linear Discriminant Analysis; 
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Introduction 

Edible oils bring essential nutrient 
components for human beings such as vitamins, 
fatty acids, and micronutrients, which are 
necessary for daily life [1]. However, the use of 
expired edible oils leads to a decrease in the 
nutritive value and an increase in potential 
hazards to people's health [2]. Therefore, the 
authentication and identification of edible oils 
are of great importance in the field of food 
safety and quality monitoring. 

A reliable, fast and non-destructive detection 
method to identify various types of edible oils 
is essential. Recently, different methods have 
been employed to identify the types of edible 
oils [3]. Several works have been reported in 
literature exploring the use of liquid 
chromatography [4], fluorescence spectroscopy 
[5], fourier transform infrared [6], differential 
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scanning calorimeter [7], supercritical liquid 
chromatography (SFC) and gas 
chromatography [8], for the quality control of 
vegetable oils. However, most of these methods 
are highly sophisticated, expensive and involve 
laborious analysis. 

Lipid oxidation in vegetable oils is 
associated with unsaturation of the oils. This 
reaction leads to the formation of a series of 
intermediate compounds named hydro 
peroxides. Hydro peroxides are the primary 
oxidation products of lipid oxidation. This is 
due to the unstable nature of these products 
which leads to their decomposition and turning 
into secondary oxidation products, such as 
carbonyl compounds. These products are 
generally unstable and decompose into a 
variety of secondary oxidation products, 
including carbonyl compounds [9]. 
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This condition causes flavor and odor 
change; Hydro peroxides do not play an 
important role in flavor deterioration, whereas 
carbonyl compounds are mainly responsible for 
the typical rancid off-flavors [10]. 

Two of the best analytical indicators of 
oxidative changes in fats and oils are peroxide 
(PV) and Carbonyl (CV) values [11]. 
Determination of these values provides 
valuable information in regard to primary and 
secondary oxidations [12]. However, due to the 
use of chemical materials (solvents such as 
chloroform, n- Hexane etc.), determination of 
these values is time- consuming, financially 
costly and costly in terms of damage to 
environment. Therefore, finding a quicker, less 
expensive and more environmentally friendly 
method is crucial. 

Nowadays, digital imaging is becoming 
more important because of its ability to perform 
fast and non- invasive, low- cost analysis on 
foods. In fact, a wide variety of digital cameras 
and digitalization equipment has contributed to 
an increase in the number of papers published 
exploring the use of webcam [13], scanner 
[14], cell phones [15] and digital camera [16] to 
monitor the quality of several food samples. A 
positive feature of using digital image to 
monitor the quality of foods is that it replaces 
the human visual system, often employed in 
these types of analyses [17]. 

Therefore, the use of digital image 
eliminates the subjective character of analyses 
as well as the dependence of the human visual 
system, which is substantially influenced by 
ambient conditions and _— subject to 
inconsistencies [18]. 

Given this, techniques based on digital 
images are a promising alternative for the 
analysis of food and other products. Some 
studies have shown the availability of using the 
technique for quality control in various 
matrices, including shrimp [19], cereal grains 
[20], Kiwifruit [21], castor seeds [22], olive oils 
[23-25] and other edible vegetable oils [26]. 
Most of these studies use information from 
color models associated with pattern 
recognition techniques to cluster, or classifythe 


samples into categories, according to similarity 
standards. 

The red, green and blue additive color added 
together to form RGB color space that was 
designed to match an intuitive human 
perception of the colors. In HSI color space, 
parameters are the intensity I, chromaticity hue 
H and saturation S. HSI approximates the way 
in which humans perceive and interpret colors. 

At present, usually, the color of foods has 
been measured in L*a*b*. The L*a*b*, or CIE 
L*a*b*, color space is an international standard 
for color measurements, adopted by the 
Commission Internationale d’Eclairage (CIE) 
in 1976 [27]. 

However, at present available commercial 
colorimeters, measure L*a*b* only over a very 
few square centimeters, and _ thus their 
measurements are not very representative in 
heterogeneous materials such as most food 
items [27]. Some of the instruments most 
frequently used in the measurement of color are 
colorimeters (e.g Minolta chroma meter, 
Hunter Lab colorimeter and Dr. Lange 
colorimeters). There is, however, a 
disadvantage in using them [26, 28-29]. 

Using a computational technique with a 
combination of a digital camera, and image 
processing software has already been used to 
provide a less costly and more adaptable way to 
measure the color of many food products and 
foodstuffs instead of traditional color 
measuring instruments [27]. 

Although there have been many studies in 
the field of image processing and classification 
of foods, studies on the classification of edible 
vegetable oils based on image processing have 
been very limited especially in the field of 
determination of expiration status (rejected and 
non-rejected). 

As was pointed before, understanding the 
first oxidation based on odor and flavor is not 
possible since the off flavor occurs in the 
second oxidation, therefore, it is critical to find 
a new technique for discerning oxidation. This 
is possible based on the color of oil, however, 
discerning the condition of rancid is impossible 
by the naked eye, therefore, it is necessary to 
use pattern recognition and image processing. 
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According to no significant color change just 
before and after rejection, the aim of this study 
was to determine rejected and non-rejected 
status of edible vegetable oils during oxidation 
process in specified intervals of two and six 
hours according to the induce period of 
peroxide and carbonyl value by using only the 
oil color changes instead of expensive and 
laborious current recognition and quality 
control methods of oils. To the best of our 
knowledge, this is the first study to investigate 
the expiration status of edible vegetable oils 
according to exact rejection point by using the 
color changes. 


Materials and Methods 

Four vegetable oil samples _ including 
Soybean, Sunflower, Palm, and Olive oils were 
purchased from the Segol Co. (Nishaboor, 
Iran). They were stored at -18°C until the time 
of analysis. All chemicals and solvents were 
purchased from Sigma-Aldrich (St. Louis, MO) 
and Merck (Darmstadt, Germany). Images were 
taken in every two (Soybean, Sunflower oils) 
and six hours (Palm, and Olive oils), during the 
heating process (at 85°C in an oven). Also, the 
chemical (PV and CV) analysis was performed 
at the same time of imaging. 


Peroxide value 

The spectrophotometric method of Shantha 
and Decker [30] was used to determine the 
peroxide value (PV). The oil samples were 
mixed with chloroform: methanol in a glass 
tube. Ammonium thiocyanate and iron (ID 
chloride was then added. After that the sample 
was placed for 5-min at room temperature, the 
absorbance of the sample was read at 500 nm 
against a blank by using a spectrophotometer. 
Results were expressed in milliequivalents of 
oxygen per kilogram of oil. 


Carbonyl value 

The carbonyl value (CV) was measured by 
using 2-propanol and 2, 4-decadienal as solvent 
and standard, respectively. Any trace of 
carbonyl compounds which may present in the 
solvent was removed with mixed solution of 


sodium borohydride and 2-propanol. Standard 
aldehyde (2, 4-decadienal) dissolved in the 
solvent. The oil was mixed with 2, 4- 
Dinitrophenylhydrazine (DNPH) solution in a 
test tube. Then the test tube was stoppered and 
heated. After that, the sample was cooled in 
water, and 2% KOH solution was added. After 
centrifugation, the absorbance of the upper 
layer was measured at 420 nm by using a 
spectrophotometer against a blank. The results 
were expressed in micromols of 2, 4 decadienal 
per gram of oil [31]. 


Apparatus and digital image acquisition 

Apparatus 

The size of system that was built for imaging 
was 120 cmx 90 cmx 90 cm with dark walls to 
isolate the samples from external light. The 
compartment has camera canon model, EOS 
1000D, which was connected to computer by a 
USB port. The illumination of the compartment 
was performed by using eight fluorescent lamps 
with 8 W (white color), the lamps were placed 
at a distance of 20 cm and 40 cm from the 
sample and with 45 degree angle for preventing 
the reflection of light. 

The [Illustration was performed by 
Zoombrower EX 0.5, the other characteristics 
of camera for imaging were as follow: flash 
(off), zoom (on), Iso speed (100), Aperture 
priority (F / 20) and Shutter speed (0.6 Sec). 
The illumination condition at compartment for 
each sample was the same. 


Image color analysis 

Image color analysis was performed using 
the MATLAB (R2013) software to convert 
images from R*G*B color space to L*a*b, HSI, 
Grayscale. 

The recorded images contained 24-bit (16.7 
million colors) and 3888 pixelsx 2592 pixels 
spatial resolution, were stored in JPEG format 
QGpg). A circular region with a radius of 2.75 
mm at the center of each image was selected for 
converting R*G*B to L*a*b, HSI and 
Grayscale. The extracted color values were then 
used for linear discriminant analysis (LDA) 
classification (Fig. 1). 
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Fig. 1. System used in the acquisition of images and classification of oils samples. 


Chemometric procedures 

The extracted color features were analyzed 
based on two approaches. The first approach 
involved the classification of non-rejected and 
rejected status of each oil samples with respect 
to induced period of peroxide value, and the 
second classification involved the classification 
with respect to induced period of carbonyl 
value, to distinguish non-rejected and rejected 
status of edible oils according to primary and 


secondary oxidation. A total of 99 images of 
oils data were collected during the heating 
process. Data analysis was performed by 
classical multivariate procedures including 
LDA with MATLAB (R2013). The RGB, HSI, 
L*a*b* and Grayscale color space features 
were extracted from images of oils samples 
then they were used in LDA classification for 
both color changes during primary (90 images) 
and secondary (99 images) oxidation (Table 1). 


Table 1. The total number of oil samples in both PV and CV tests. 


Palm Olive Soybean Sunflower 
Non-rejected 20 16 7 7 
PV Rejected 5 10 12 13 
Total 25 26 19 20 
Non-rejected 24 20 11 11 
CV Rejected 10 6 8 9 
Total 34 26 19 20 
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To find the effect of oxidation process on the 
color of edible vegetable oils during time, the 
color changes in each oil type at each color 
spaces RGB, HSI, L*a*b* and Grayscale was 
considered. 

Linear Discriminant Analysis (LDA) is a 
classification method that is commonly used as 
dimensionality reduction technique for pattern- 
classification and machine learning to find a 
linear combination of features that 
characterizes or separates two or more classes 
of objects which was developed by Fisher [32]. 
Multi-class LDA is a generalization of standard 
two-class LDA that can handle arbitrary 
number of classes. We are seeking p projections 
(1.e.Y1, Y2... Yp) of the input vector x by means 
of p projection vectors w; as follows: 

y, = wx, be ae CE) 


If we arrange all w; in a projection matrix as 
W=[W1,W2,...,Wp], it can be written: 


y=W'x (2) 
where 

x=[% Xm) 

and 

y=([V¥° Mni. 


By stacking all feature vectors in one matrix, 
we can write: 


Y=Ww'xX (3) 
where 
X = [x1,x?, a| a | 
Ming Ke: es Be 
(4) 
Mw v. M 
Valparaiso. SA 
Va Ve Dy 
(5) 


For p-classes case, we will measure the 
within-class and between-class scatters with 
respect to the mean of all classes, respectively, 
as follows: 


Sw = iat eee (k= By) ep, )? (6) 
Sp = Die Nilmi — Bw) (Hi — BIT (7) 


Where 
1 1up 
M= dv, X = Fda Mimi (8) 
1 
Bj = nrc © (9) 


Such that C; represents the i-th class. 
Similarly, we can define the mean vectors for 
the projected samples as follows: 


~ 1 


hi =—d y (10) 


Nj, yYECi 


Bt 

w=2y,, y (1) 

Similarly, the within-class and between- 
class scatter matrices of the projected samples 
can be given, respectively, by: 


Sw = Lie1Lyec(y — Bly — H)? (12) 
Sp = Deer Nil; — Mi) (a; — a)? (13) 


By using Eqns. (6) and (7) and some 
algebraic manipulations, we can obtain: 


Sw = W'SyW (14) 
Sp = W'S,W (15) 


To obtaine an appropriate discrimination 
between all classes, the coefficients of W should 
be optimally adjusted such that the between- 
class scatter increases while simultaneously, 
the within-class scatter decreases. Thus, the 
optimal W’ can be obtained by maximizing the 
following objective function: 


jw) = Mal. = ws 


Sal 
lSwl |wTsww| 


(16) 
W*= index(max J) (17) 


Where]|.| computes the determinant of a 
matrix. For W’, it is sufficient to set the 
differential of J (W) (with respect to W) equal 
to zero. It can be shown that the columns of W” 
are the eigenvectors corresponding to the 
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largest eigenvalues of the following generalized 
eigenvalue problem: 


Sy Spw; = A,w; (18) 
Where 


Wi = [Wt, W2,...Wp-1]. 
Evaluation Measure 
This feature was studied by using the 
accuracy (ACR) analysis as follow: 
eps TP+IN (19) 
TP+FP+TN+ FN 
Where TP, TN, FP, and FN are the true 


positive, true negative, false positive and false 


negative, respectively, TP (FP) or true positive 
(false positive) means all samples correctly 
(incorrectly) identified. Similarly FN or false 
negative means that those samples are 
incorrectly identified. 


Results and Discussion 

Table 2 illustrates the Chemical indicators, 
Peroxide values (PV) and Carbonyl values 
(CV) of the oil samples. The PV and CV values 
of the soybean, sunflower, olive and palm oils 
after heating process are shown in Table 2. The 
PV and CV of the non-rejected oils were lower 
than 2 meq O2/kg oil and 3.1 umol/g for all 
samples and their levels after heating did differ 
significantly to reach their rejection points. 


Table 2. Chemical indicates and Peroxide value (PV), and carbonyl value (CV) of the oils before and 
after heating process at 85°C. 


Pv? IP,vTime (hour) CV‘ IP.,Time (hour) 

Palm R 0.69+40.021 3.06+40.012 
Olive R 1.99+0.46 3.047+40.07 

2 
Sunflower R 0.28+40.034 0.26+40.005 
Soybean R 0.82+0.057 0.136+40.00 

1 
Palm NR 62.74+0.11 119.09 8.74+40.21 140.68 
Olive NR 32.05340.90 93.64 10.9640.15 114.66 
Sunflower NR = 2.23+0.13 13.70 4.61+40.03 21.42 
Soybean N 1.8940.19 12.94 4.52+0.04 21.39 


a Mean value + standard deviation, All values are means of three determinations. 


b Peroxide value (meq O2/kg oil). 
c Carbonyl value (umol/g). 


During the heating process, a wide range of 
PV and CV was observed among the oil 
samples every 2 and 6 hour (from 0.82 to 1.89 
meq O2/kg oil for Soybean, 0.28 to 2.23 
meq O2/kg oil for sunflower, 1.99 to 32.053 
meq O2/kg oil for olive and from 0.69 to 62.74 
meq O2/kg oil for palm and also from 0.136 to 
4.52 umol/ g for Soybean, 0.26 to 4.61 umol/ g 
for sunflower, 3.047 to 10.96 umol/ g for olive 
and from 3.06 to 8.74 umol/ g for palm, PV and 
CV values ,respectively), indicating that the 
oils had no similar rejection points. Table 2 also 


demonstrates the IPPV and IPCV for each oil 
type. The peroxide value (PV) for soybean was 
12.94, sunflower 13.70, olive 93.64 and palm 
119.09 hour And similarly carbonyl value (CV) 
for soybean was 21.39 sunflower 21.42, olive 
114.66 and palm 140.68 hour. 

These chemical analyses were carried out in 
order to find the exact rejection point of primary 
and secondary oxidation of these oils, so that 
their classification could be done accurately in 
all their different rancidities (Fig.2). 
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Fig. 2. A schematic kinetic curve of peroxide and carbonyl accumulation during oxidation of lipid systems in 
Sunflower and Olive oil. 


By increasing the amount of PV and CV suggested that the oxidation can significantly 
during time of heating, the color of the oils affect the color. 
changed in each stage. According to the In classification part, Table 3 illustrates the 
difference in values and color feature extraction classification accuracy of non-rejected and 
of images, it may be possible to make a model rejected of each known oil type during the 
for classification of expiration status of these heating time at each color space separately to 
oils without the necessity to perform high cost identify the expiration status according to 
chemical experiments. Therefore, it is carbonyl and peroxide values. 
Table 3. Accuracy of each known oil type during the heating time at each color space in both PV and 
CV tests. 
Olive Palm Soybean Sunflower 
L*a*b* 84.61% 100% 89.48% 90% 
PV RGB 84.61% 96% 78.95% 100% 
HSI 80.77% 100% 84.21% 100% 
Grayscale 80.77% 96% 95% 715% 
L*a*b* 100% 100% 84.21% 90% 
cv RGB 96.15% 91.18% 78.95% 90% 
HSI 84.62% 88.25% 89.47% 95% 


Grayscale 88.46% 100% 78.95% 10% 
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As it can be seen in Table 3, all rejected and 
non-rejected oils according to peroxide induced 
period in each color space and Grayscale were 
identified appropriately. It is apparent from the 
Table, that the highest accuracy among studied 
color spaces were 100% at L*a*b* and HSI for 
palm, 84.61% at L*a*b*, RGB color spaces for 
olive, 95% at Grayscale color space for soybean 
and finally 100% at RGB and HSI color spaces 
for sunflower. 

Similarly, Table 3 shows the accuracy of 
edible oils according the carbonyl induced 
period. This indicated that, the best results for 
palm, olive, soybean and sunflower oils were 
achieved with 100% at L*a*b* and Grayscale, 


100% at L*a*b*, 89.47% at HSI and lastly 95% 
at HSI color space respectively. 

Table 4 illustrated the classification 
accuracy of four oils (Palm, Olive, Soybean and 
Sunflower) in non-rejected and rejected status 
in both oxidations. At first all non-rejected oils 
were classified in each color space (L*a*b*, 
RGB, HIS and Grayscale). As shown in Table 
4 the best result belongs to HSI with 88% in 
primary oxidation and with regard to secondary 
oxidation L*a*b* with 91.04% have the best 
accuracy. For rejected oils classification the 
best color space are HSI and L*a*b* with 90% 
and 87.5% in primary and secondary oxidation, 
respectively. 


Table 4. Accuracy of Non-rejected and Rejected four oils (Palm, Olive, Soybean and Sunflower) 
during the heating time at each color space in both PV and CV tests. 


L*a*b* RGB HSI Grayscale 
PV Non-rejected oils 86% 80% 88% 68% 
Rejected oils 80% 77.5% 90% 715% 
cv Non-rejected oils 91.04% 68.66% 82.09% 58.21% 
Rejected oils 87.5% 81.25% 71.88% 56.25% 


Figs. 3A and 3B show the score plots of the 
two discriminant functions (F2 * Fl) obtained 
by the LDA classifier for the R and NR oils 
samples. It is possible to observe that Soybean- 
R, Sunflower-R, Olive-R and Palm-R oil 
samples are separated along the F1 direction. It 
is obvious that the Fl separated Soy-R and 
Palm-R groups from other groups. 

As can be seen, the Olive-R and Sunflower- 
R sample has overlap with each other (Fig. 4A). 
Fig. 4B illustrated that Sunflower-NR, 
Soybean-NR, Olive-NR and Palm-NR oil 
samples presented a tendency of separation 
along the F1 direction. Although Sunflower- 
NR and Olive-NR samples well discriminated 
but Palm-NR and Soybean-NR have a slight 
overlap. 

As it can be seen in Table 5, all rejected and 
non-rejected oils according to both oxidations 
in all color spaces and Grayscale were 
classified. It is apparent from Table 5 that the 
highest accuracy belongs to the classes with 


more sample numbers (according to Table 2). 
For example, Olive NR with 100% in both 
oxidations has best accuracy in comparison 
with the Olive R (30% and 50%). Because of 
different rate of expiration in each oil-type, the 
number of samples is quite limited in some 
classes. However, total accuracy of LDA 
classifier in primary and secondary oxidation 
relatively well performs with 74.44% and 
75.16%, respectively. 

The classification performance of LDA 
classifier is demonstrated by Fig. 4A and B, 
which presents the score plots of the first two 
discriminant functions (F2 * Fl) for all Non- 
rejected and Rejected oils during the heating 
time at all color spaces in each oxidation. In 
case of carbonyl test (Fig. 4a), it is obvious 
Olive oil samples in both status (NR and R) are 
separated along the F2 direction with slightly 
overlap. Fl clearly distinguishes Olive oil 
samples from other oils classes. 


Discerning expiration status of edible vegetable oils based on color changes ... 153 


| # 
fei 2 au 
x fe A A 
2 | & | = -. 
= e T re 4 T 
a 5 
= Se 
23 + of 
e 
-8 
-9 -4 1 6 
F1 (88.14 %) 
7 
B 
5 m2 
3 e x 


F2 (5.37 %) 
& 
+ 
on) 
4 
c 
» 
®o,. 
| ae 
& 
Sg 
aa 
id 


=l a Ag . 
, ca 
ag ¢ ae 
Est 
“§ of 
7 
-7 “5 3 A l 3 5 7 

F1 (94.53 %) 


Fig. 3. F2x F1 score plots for Non-rejected and Rejected oils during the heating time atL*a*b*color space in 
carbonyl test (a) All R oilsand (b) All NR oils. Olive; @: Palm; A: Soybean; Mand Sunflower; © 
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Fig.4. F2x F1 score plots for all Non-rejected and Rejected oils (Palm, Olive, Soybean and Sunflower) during 
the heating time at all color spaces (a) carbonyl test and (b) peroxide test. Olive NR; ©: Olive R; 0: Palm NR; A: 
Palm R; A: Soybean NR; ": Soybean R; O: Sunflower NR; @and Sunflower R; ©. 
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In the case of Sunflower oil samples, this 
class was well discriminated, but the NR and R 
status are intense overlap in both directions. 
The Palm oil samples overlap with other oil 
classes but the NR and R status were separated 
well in F1 direction. Finally, the worst 
classification performance was presented by 
Soybean oil samples in both directions which 
had severe overlap with Palm-NR. 

Milanez and Pontes [25] tried to distinguish 
the type and conservation state of four different 
edible oils based on color image processing and 
linear discriminant analysis (LDA). They only 
attempted in classification of different type of 
oil, and also classification expired and non- 
expired oils without considering the oxidation 


process of oils during time with analytical 
experiments. 

It could be inferred from our study that, it is 
appropriate to use color space to identify the 
expiration status of each oil type during heating 
process. As muchas the oil color changes with 
different rates during oxidation, it could be 
possible to detect the exact rejection state 
(based on both the primary and secondary 
oxidation) of known oil sample. 

It seems that while classification of edible 
oils from color features can be applicable, 
selection of the suitable color space and 
classifier are significant steps to develop new 
simple and inexpensive method for detection of 
the non-rejected and rejected status of unknown 
oil. 


Table 5. Accuracy (%) of all Non-rejected and Rejected oils (Palm, Olive, Soybean and Sunflower) 
during the heating time at all color spaces in both PV and CV tests. 


Olive Olive Palm Palm 


Soybean Soybean Sunflower Sunflower Sample 


NR R NR eR NR R NR R Number ccuracy 
PV 100 30 93.33 90 63.64 62.50 44.44 81.82 90 74.44 
CV 100 50 92.31 37.50 40 44.44 81.82 88.89 99 75.76 
Conclusions Grayscale), olive with 100% (L*a*b*), soybean 


The purpose of this paper was to find less 
costly and quicker method with environmental 
protection, by using the color spaces (RGB, 
HSI, L*a*b* with Grayscale) instead of 
chemical analysis to determine the expiration 
status of edible vegetable oils in both PV and 
CV induced period. 

The results obtained from this study 
indicated that the best classification result of 
expiration status of known oils according to 
induced period of peroxide value at each color 
space, was achieved with LDA model were in 
palm with 100% (L*a*b* and HSD), olive with 
84.61% (L*a*b* and RGB), soybean with 95% 
(Grayscale) and sunflower with 100% (RGB 
and HSD), also in induced period of Carbonyl 
value test, the best classification performance 
was achieved in palm with 100% (L*a*b* and 
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